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ABSTRACT 

We present our NextGen Model Atmosphere grid for low mass stars for 
effective temperatures larger than 3000 K. These LTE models are calculated 
with the same basic model assumptions and input physics as the VLMS part 
of the NextGen grid so that the complete grid can be used, e.g., for consistent 
stellar evolution calculations and for internally consistent analysis of cool star 
spectra. This grid is also the starting point for a large grid of detailed NLTE 
model atmospheres for dwarfs and giants (Hauschildt et al, in preparation). 
The models were calculated from 3000 K to 10000 K (in steps of 200 K) for 
3.5 < \og{g) < 5.5 (in steps of 0.5) and metalhcities of -4.0 < [M/H] < 0.0. 

We discuss the results of the model calculations and compare our results to 
the Kurucz 1994 grid. Some comparisons to standard stars like Vega and the 
Sun are presented and compared with detailed NLTE calculations. 
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1. Introduction 

As part of our ongoing effort to calculate accurate line blanketed stellar atmospheres we 
present the results our "NextGen" model atmosphere grid for "high" effective temperatures 
(3000 < Teff < 10000 K). The results for very low mass stars (VLMS) and Brown 
Dwarfs involve the complication of dust condensation and are discussed in AUard et al, 
(in preparation). All the NextGen models are computed with the same assumptions to 
allow internally consistent analysis of, e.g., HR diagrams of globular clusters or population 
synthesis models. 

Most of the models have been calculated under the assumption of local thermodynamic 
equilibrium (LTE) for the ionization/dissociation equilibria and the level populations. 
Thomson and Rayleigh scattering are included and cause the source function to deviate 
from the Planck function. The LTE assumption has been made mainly to be consistent with 
the VLMS part of the NextGen grid. The LTE models are being used as starting points 
for a large NLTE model grid that is currently being constructed and that will replace the 
LTE model grid (Hauschildt et al, in preparation). The NLTE models that are currently 
available in the NextGen model grid feature only a small subset of NLTE species and will 
be replaced by full NLTE models. We will discuss briefly how NLTE models change the 
results for one example. 

Kurucz (1992) has calculated an extensive grid of LTE model atmospheres with his 
ATLAS9 code using Opacity Distribution Functions (ODE) to include line blanketing. This 
model grid has set the standard for a large number spectral analyses over a wide range of 
effective temperatures. Castelli & Kurucz (1994) have used the Kurucz 92 grid and the new 
opacity sampling version of Kurucz's model atmosphere code ATLAS12 (Kurucz, 1993a) to 
compute model atmospheres for Vega. They find very good agreement over a wide range in 
wavelengths, improving significantly upon earlier results. This is mostly due to Kurucz's 
update to his line lists, resulting in an increase in the total line opacity resulting in better 
fits to the ultraviolet spectra. 

We use the Kurucz 92 grid to test our calculations at effective temperatures larger than 
about 5000 K. At lower temperatures, larger differences are expected due to differences 
between the equations of state (EOS) and the treatment of molecular opacities in ATLAS 12 
and PHOENIX. For higher effective temperatures, when molecules and molecular opacities 
become less important, the synthetic spectra should be in good agreement with our LTE 
models because we use the Kurucz atomic fine list (Kurucz, 1994); small differences are 
expected due to different bound-free and free-free opacity sources and different numerical 
methods used. 
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2. Methods and Models 

For the model grid calculations, we use our multi-purpose stellar atmosphere code 
PHOENIX (version 9.1, Hauschildt et al., 1996; Hauschildt, Baron, & Allard, 1997; Hauschildt 
et al., 1997a; Baron & Hauschildt, 1998). In the calculations presented in this paper, we will 
mostly use the static, plane-parallel LTE mode of PHOENIX, although we will compare LTE 
results with NLTE calculations to estimate the systematic errors made by LTE analyses. 
Details of the numerical methods are given in the above references, so we do not repeat 
the description here. We will, however, describe the changes in the input physics and data 
compared to our earlier "base" and "Extended" model grids for very low mass stars. 

Both atomic and molecular lines are treated with a direct opacity sampling method. 
We do not use pre-computed opacity sampling tables, but instead dynamically select the 
relevant LTE background lines from master line lists at the beginning of each iteration for 
every model and sum the contribution of every line within a search window to compute 
the total line opacity at arbitrary wavelength points. The latter feature is crucial in NLTE 
calculations in which the wavelength grid is both irregular and variable from iteration to 
iteration due to changes in the physical conditions. This approach also allows detailed and 
depth dependent line profiles to be used during the iterations. 

Although the direct line treatment seems at first glance computationally prohibitive, it 
can lead to more accurate models. This is due to the fact that the line forming regions in 
cool dwarfs span a huge range in pressures and temperatures so that the line wings form in 
very different layers than the line cores. Therefore, the physics of the line formation is best 
modeled by an approach that treats the variation of the line profile and the level excitation 
as accurately as possible. In addition, it is relatively straightforward to write an efficient 
computer code for the direct opacity sampling treatment by making use of data-locality, 
vectorization & super-scalar execution, and caches found in modern workstations and 
supercomputers. On many machines, interpolation in the large tables required by the 
ODF method takes more time per wavelength point than direct opacity sampling because 
of the longer memory access times compared to the execution speed of floating point 
instructions. To make this method computationally more efficient, we employ modern 
numerical techniques, e.g., vectorized and parallelized block algorithms with high data 
locality (Hauschildt, Baron, & Allard, 1997), and we use high-end workstations or parallel 
supercomputers for the model calculations. 

In the calculations presented in this paper, we have have included a constant statistical 
velocity field, ^ = 2kms~^, which is treated like a microturbulence. The choice of fines 
is dictated by whether they are stronger than a threshold T = xi/K'c — 10~^, where xi is 
the extinction coefficient of the line at the line center and Kc is the local b-f absorption 
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coefRcient. This typically leads to about 10 x 10^ lines which are selected from master hne 
lists (with 47 million atomic and up to 350 million molecular lines). The profiles of these 
lines are assumed to be depth-dependent Voigt or Doppler profiles (for very weak lines). 
Details of the computation of the damping constants and the line profiles are given in 
Schweitzer et al. (1996). We have verified in test calculations that the details of the line 
profiles and the threshold F do not have a significant effect on either the model structure or 
the synthetic spectra. In addition, we include abound 2000 photo-ionization cross sections 
for atoms and ions (Mathisen, 1984; Verner & Yakovlev, 1995). 

2.1. Moleculcir line data 

One of the most important improvements of the new models compared to the Allard & 
Hauschildt (1995) models (hereafter, AH95) is that new molecular line data have recently 
become available. The first improvement of our molecular hne list is the addition of the 
HITRAN92 (Rothman et al., 1992) database, first incorporated in our "Extended" model 
grid (Allard & Hauschildt, unpublished). For the molecules with lines available from other 
sources (e.g., water vapor and CO) we prefer to use other sources in model calculations 
because the HITRAN92 lists are fairly incomplete (they have been prepared for the conditions 
of the Earth atmosphere) although the line data are of very high quality. The OH molecule 
is an exception because the vibrational bands present in the HITRAN92 database are not 
present in the Kurucz CD 15 (Kurucz, 1993b) datasct listing electronic transitions of OH. 
Therefore, we use both HITRAN92 and CD 15 OH line data simultaneously in the model 
calculations. 

We have replaced the vibration-rotation lines of CO (Kurucz, 1993b) used in AH95 with 
the more recent CO line list calculated by Goorvich (Goorvitch & Chackerian, 1994a,b). 
This hst is more accurate than the old CO hne data. However, we have kept the electronic 
CO transitions available on Kurucz's CD 15. Similarly, we have replaced the Kurucz CN 
line list (CD 15) with a more recent calculation by J0rgensen (J0rgensen & Larsson, 1990). 
CN lines are comparatively weak in most of the models presented here. In addition, we have 
added line lists for YO (Littleton, 1987) and ZrO (Littleton & Davis, 1985) as well as H^ 
partition functions and lines (Neale & Tennyson, 1995; Neale, Miller, & Tennyson, 1996). 

Another improvement over the "Base" model grid (Allard & Hauschildt, 1995) is the 
replacement of the straight mean TiO opacities with the list of TiO lines computed by 
J0rgensen (J0rgensen, 1994). This list includes a total of more than 12 million lines for the 
most important band systems of 5 TiO isotopes. On average, the opacity of TiO lines is 
about a factor of 2 smaller than the old straight mean opacities. For any given effective 
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temperature this will result in weaker TiO features when compared to the old straight 
means. 

The absorption coefficient of water has been a long standing problem for modelers of 
cool dwarf atmospheres. In our Base model grid we have used the straight means of the 
Ludwig (1971) water opacity tables. These are known to overestimate the water opacity 
significantly at higher (and for M dwarfs more important) temperatures (Schryber, Miller, 
& Tennyson, 1995) but they are accurate at lower temperatures. The currently available 
list of water lines of Miller & Tennyson (MT-H2O, Allard et al., 1994) is incomplete and 
therefore some water opacity will be missing in models with Tcs < 4000 K. Comparison with 
high-resolution observations in spectral regions where the MT-H2O list should be complete 
show that its accuracy is very high (Jones et al., 1995). In the models presented here, 
we have used the MT-H2O fist. Recently, a new water vapor fine list with more than 300 
million lines was published by Partridge & Schwenke (1997). A new and improved version 
of the Miller & Tennyson list is currently being calculated (Viti et al, in preparation). In a 
subsequent publication wc will show the differences in the model structure and the synthetic 
spectra for low effective temperature models calculated using different water vapor line lists. 

2.2. The equation of state 

The equation of state (EOS) is an enlarged and enhanced version of the EOS used in 
AH95. Wc include about 500 species (atoms, ions and molecules) in the EOS. This set of 
EOS species was determined in test calculations. The EOS calculations themselves follow 
the method discussed in AH95. We did not include dust condensation in the calculations 
presented here because it does not affect the models in the range of parameters we consider. 
Models for effective temperatures lower than about 2500-3000 K need to include these 
effects and will be discussed in a subsequent paper. Some of the molecular data were 
updated, details are given in the aforementioned paper. 

3. Results 
3.1. The Grid 

The NextGen model grid was calculated from 3000 K to 10000 K (in steps of 100-200 K) 
for 3.5 < log(^) < 5.5 (in steps of 0.5) and metaUicities of -4.0 < [M/H] < 0.0 (in 
steps of 0.5). We use the recent solar abundances given in Table 5 of Jaschek & Jaschek 
(1995). The changes in the abundances compared to the previous Anders & Grevesse 
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(1989) data are in general small, with the exception of the iron abundance which changed 
from 7.67 to 7.53. The grid contains a total of 2142 models in this range of parameters 
(the complete NextGen grid includes 3026 models). In all models, convection is treated in 
the mixing length approximation with the mixing length set to unity. We will discuss the 
model grid for lower effective temperatures, the VLMS to Brown Dwarf range of spectral 
class M, in a separate paper. For effective temperatures higher than 7000-10000 K NLTE 
effects become important, so LTE models at higher effective temperature will become 
increasingly unrealistic. In cooler models, NLTE effects are important for individual atomic 
lines of, e.g., Ti I (Hauschildt et al., 1997b), for effective temperatures below about 
4500 K but NLTE does not have significant effects on the low resolution spectra or the 
structure of these models. Therefore, we include a few representative NLTE models (for 
7000 < Tcflf < 10000 K) in the model grid presented here. A full NLTE grid (with the 
complete set of NLTE species available in PHOENIX) is currently being constructed and will 
be discussed in a later paper (Hauschildt et al, in preparation). 

Aufdenberg et al. (1998a,b) have shown that for effective temperatures higher than 
about 18000 K the combined effects of line blanketing and spherical geometry are of crucial 
importance for main sequence stars with log(5f) < 4.5. In this parameter range, plane 
parallel models do not deliver enough EUV flux compared to spherical models or observed 
EUV spectra. However, for effective temperatures below about 10000 K the use of the 
plane parallel approximation in the model construction for log(g') > 3.5 does not result 
in significant differences to spherical models. Therefore, we have used the plane parallel 
geometry to calculate the present grid. A grid of NLTE spherical model atmospheres for OB 
main sequence stars is currently being constructed and will be discussed in a separate paper 
(Aufdenberg et al, in preparation). For effective temperatures lower than about 3000 K the 
effects of dust formation and/or dust opacity become important. This significantly changes 
the physics of the model atmospheres and the formation of the spectrum, therefore, we will 
discuss these models in a separate paper (Allard et al, in preparation). 

The synthetic spectra and the model structures are available via anonymous 
FTP from ftp://calvin.physast.uga.edu/pub/NextGen or via the WWW URL 
http://dilbert.physast.uga.edu/~yeti and constitute about 560 MB of data. For 
each model, the model structure (in the form of the PHOENIX output file of the last 
model iteration) and the synthetic low resolution spectra (directly the result of the model 
iterations) are available. The model fluxes are given as tables of Fx in erg/s/cm^/cm versus 
wavelength in A to make comparisons with observed spectra and other model calculations 
easier. The spectra are given on an semi-regular wavelength grid with about 23,000 
wavelength points from lA to 1 cm, the actual wavelength grid depends on the model. 
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3.2. Compcirison to Kurucz 92 models 

In Figs. 1-3 we compare the PHOENIX synthetic spectra with ATLAS9 spectra for a range 
of effective temperatures with fixed solar abundances. Note that while Kurucz 92 uses the 
Anders & Grevesse (1989) solar abundances, we use the more recent solar abundances given 
in Table 5 of Jaschek & Jaschek (1995). The changes in the abundances will introduce 
small but systematic differences between the Kurucz 92 models and our model structures 
and between the spectra because the abundances of the key elements C, N, O, and Fe are 
different. The PHOENIX spectra have been degraded by convolution with a Gaussian kernel 
to more closely resemble the resolution of the ATLAS9 spectra. The spectra for the Kurucz 
models were obtained using the IDL routine KURGETl and the corresponding database of 
models available in the lUE reduction and data analysis package lUERDAF. 

In general, our models agree well with the ATLAS9 spectra for effective temperatures 
higher than about 5000 K. We expect better agreement between PHOENIX and ATLAS9 
spectra for higher effective temperatures since we have included the latest version of 
Kurucz atomic line list into PHOENIX and atomic lines dominate the overall opacity in this 
temperature range. The remaining differences between the spectra are probably caused by 
the shghtly different abundances, different methods for treating the line blanketing (ODF in 
ATLAS9 versus direct opacity samphng in PHOENIX), the fact that the ODF's in the ATLAS9 
models were constructed using a previous version of the Kurucz line lists, and by differences 
in the treatment of, e.g., scattering and b-f and f-f opacities. These different treatments 
alter the model structures, resulting in systematic variations between the spectra. In the 
ultraviolet (UV), the two model grid spectra differ in several regions (cf. Fig. 1), which 
might in part be due to the fact that the ATLAS9 spectra are based on ODF's. The general 
features of the UV spectra are similar. Some of the features in the PHOENIX spectrum are 
deeper than the ATLAS9 features. This could partly be due to the different abundances and 
to changes in the atomic line list itself since the calculation of the ATLAS9 models. 

We compare the models in the region of the H I Balmer lines in Fig. 2. The lines widths 
agree well, however, the H I lines in the PHOENIX spectrum seem to be a little deeper than 
the ATLAS9 lines. This is probably simply an artifact of the ODF versus opacity sampling 
treatment of the lines. In addition, we have inserted wavelength points in the cores on the 
H I lines. For comparison, we also show the PHOENIX spectrum at its original resolution with 
the fluxes scaled by a factor of 0.5. The opacity samphng method results in many more hue 
features appearing in the spectrum, similar to ATLAS 12 results shown in Castelli & Kurucz 
(1994). Figure 1 compares our models to Kurucz 92 for log(5f) = 4.0, Tcs = 10000 K and 
for solar abundances. The differences between the spectra are essentially the same over the 
whole range of gravities and are therefore not shown. The PHOENIX spectra show less flux 
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in some UV regions (e.g., around 1500A) and the UV lines are generally deeper than in the 
ATLAS9 spectra. The flux that is intercepted in the UV is redistribTitcd into the Paschen 
continuum so that the total energy flux through the atmosphere is given by aTcfi'^. Thus, 
the differences in the spectra are probably caused mainly by a larger UV opacity in the 
opacity sampling PHOENIX models. 

For lower effective temperatures, the differences between the PHOENIX and the ATLAS9 
spectra become larger. In Fig. 3 we show a more detailed comparison between models with 
Tcs = 4000 K, log{g) = 4.0 and solar abundances. The figure clearly shows the spectral 
regions where the two spectra do not agree well. In particular, the "pseudo-continuum" 
between about 0.9 ii and 1.6/i lies at about a 10% lower flux level in the ATLAS9 spectrum 
than in the PHOENIX spectrum. In the optical spectral region between 6OOOA and about 
8OOOA the PHOENIX spectrum shows a number of molecular bands (due to TiO and VO?) 
that are not seen in the ATLAS9 spectrum. These bands are more important in the VLMS 
and Brown Dwarf range of the NextGen grid. The Kurucz TiO line list results in much 
weaker TiO features and different band shapes than the J0rgcnscn TiO list that we use 
(Schweitzer, 1995). This is more important for lower effective temperatures and will cause 
more significant deviations between the PHOENIX and the ATLAS9 models. 

There are also discrepancies between the spectra in the near infrared (NIR) spectral 
region. The CO bands are somewhat different, which is probably caused by the ODF versus 
opacity sampling treatment of the line blanketing of the calculations. The CO line data 
that we use arc very similar to the Kurucz CO line list (Schweitzer, 1995). We include 
water lines in the model calculations and the spectral synthesis, which is not included in the 
ATLAS9 models. This causes the differences in the molecular absorption band around 2.5 /x, 
this feature is due to a strong water vapor band. Therefore, it is likely that the differences 
for low effective temperatures can be explained mainly by the different opacity data used in 
PHOENIX and ATLAS9. At effective temperatures as high as 4000 K, differences between the 
number of molecules and the molecular data used in the EOS's of PHOENIX and ATLAS9 are 
probably unimportant. At lower temperatures, however, the larger EOS of PHOENIX will 
likely increase the differences between the two sets of models. 

3.3. Comparison to NLTE models 

In order to investigate the importance of NLTE effects on the models presented in this 
paper, we have calculated small number of representative NLTE models with parameters 
roughly appropriate for the Sun and Vega (a Lyr). We include a large number of levels and 
lines in NLTE in order to make the model calculations as realistic as possible. In addition. 
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species that we do not include explicitly in NLTE are included in the LTE approximation 

in order to describe their effects on the model structure and the spectrum, see Hauschildt 
et al. (1997b) and Aufdenberg et al. (1998a) and references therein for details of the 
method. This allows us to determine which models require NLTE and for which models 
LTE provides an adequate description of line formation. In the outer, optically thin, regions 
of the atmosphere, NLTE effects will always become important (even in the absence of 
a chromospheric temperature rise), but for models of stellar photospheres, it is crucial 
to know if NLTE effects are important in the line forming regions. We have computed 
detailed NLTE models for two cases: for relatively cool, G-type stars we chose model 
parameters approximately appropriate for the Sun; and for hotter stars, where we chose 
model parameters appropriate for Vega. The models serve to demonstrate the effects of 
NLTE on the structure and spectra, we did not make any attempt to produce detailed fits 
to either the Sun or Vega. The models include a large number of NLTE species and use the 
same version of PHOENIX and the same input physics as the NextGen models. 

3.3.1. Solar- Type Stars 

For the solar-type model, we used Teg = 5770 K, log{g) = 4.44 and solar abundances. 
The only difference with the NextGen models is the use of PHOENIX's NLTE mode, which 
includes the effects of multi-level NLTE on the EOS, the radiative transfer (both lines 
and continua) and the structure of the atmosphere. This model is a very simple model 
for the solar spectrum, for example we have not used the best value for the mixing length 
parameter for the sun of a = 1.5 and have used only simple approximations for the damping 
constants of most spectral lines. Therefore, the model cannot reproduce the solar spectrum 
as well as a fine tuned model could, we only use it to demonstrate the effects of NLTE on 
solar type stars. 

Table 1 gives a complete fisting of all NLTE species currently available in PHOENIX. For 
the solar model, not all of these are important, therefore we have used only the following 
sub-set of species in NLTE: H I, He I-II, Li I, Na I, Ne I, Mg II, Ca II, Ti I-II, S II-III, 
Si II-III, C I-II, N I-II, O I-II, and Fe I-III. This resulted in a total of 4,143 NLTE levels 
and 49,324 primary NLTE transitions. 26,786 of the primary NLTE lines were treated with 
detailed Voigt profiles, the remaining 22,538 weak primary NLTE lines were treated with 
Gaussian profiles to save CPU time. In test calculations we verified that this does not 
change the results significantly. In addition, the NLTE models inchide 218,009 secondary 
NLTE lines as well 385,484 LTE background atomic lines and 1,566,441 LTE molecular 
lines, these lines were selected using the same criterion that we use for the LTE models. The 
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calculation was performed with a variable resolution wavelength grid with about 300,000 
points, including extra points that are inserted to resolve NLTE lines. We are currently 
adding Mg I, Si I, S I, Ca I, and H~ to our hst of NLTE species and we will produce 
improved models in cases where these species are significant. 

In Fig. 5 we show the comparison between the Kitt Peak Solar Atlas (Kurucz et al., 
1984) spectrum and the simplified solar model. The resolution of the Kitt Peak spectrum 
is extremely high (0.01 A). Therefore, we have reduced the resolution by convolving the 
spectrum with a Gaussian kernel with a width of lOA, shown as the full curve in Fig. 5. The 
same procedure was used to plot the synthetic spectrum (dotted line). In Fig. 6 we show 
the comparison at an enlarged wavelength scale. In this figure, the Kitt Peak spectrum and 
the PHOENIX spectrum were left at their original resolution (the resolution of the PHOENIX 
spectrum is variable, AA < 0.05A). The fit is in general good, although fine tuned model 
parameters would improve the comparison. 

The structure of the atmosphere is shown in Fig. 4 and in more detail in Fig. 7. 
The differences between the LTE and NLTE structures is very small. We show selected 
departure coefficients for the solar model in Fig. 8. In the line forming regions, typically 
between Tgtd — 10~^ and unity, the departure coefficients are very close to one and many 
lines thus form basically under LTE conditions. The departure coefficients drop rapidly 
at lower optical depths. In this region, the chromospheric temperature rise (which is not 
included in the model) would change the results drastically above the solar temperature 
minimum of about 4400 K. The wide range of different atoms and ions that are included 
in this model gives an overall impression of the average importance of the NLTE effects 
in the photosphere of the Sun. We expect that the other elements, not treated in NLTE 
in this model, will roughly follow the selected ions shown in Fig. 8. The NLTE effects on 
the photospheric structure and the line forming regions are small, therefore, LTE models 
can be used as a first approximation to model the photospheres of solar type stars. This is 
consistent with the results of Anderson (1989), who used a different numerical approach. 
The relatively high electron densities in the line forming region, which make the collisional 
rates more important in solar type stars compared to either hotter or cooler stars is the 
reason that LTE is a good approximation for solar type stars. 



3.3.2. Vega (a Lyr) 

The detailed analysis of the spectrum of Vega with ATLAS9 and ATLAS 12 models has 
been discussed by Castelli & Kurucz (1994). Here, we compare PHOENIX LTE and NLTE 
models computed with the parameters found for Vega by Castelli & Kurucz (1994) to the 
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observed Vega spectrum to show the effects of NLTE on the spectral energy distribution 
for stars with Vega's effective temperature. We have computed an NLTE model with 
Tcff = 9550 K, log(5') = 3.95 and the "Vega abundances" given in table 4 of Castelli & 
Kurucz (1994). We have used the following species in NLTE: H 1, He I-II, Ne I, Mg II, 
Ca II, Ti I-II, S II-III, Si II-III, C I-II, N I-II, O I-II, and Fe I-III, which includes the 
most important species for the conditions in Vega. This resulted in a total of 3,787 NLTE 
levels and 47,593 primary NLTE transitions. 25,449 of the primary NLTE hues were treated 
with detailed Voigt profiles, the remaining 22,144 weak primary NLTE lines were treated 
with Gaussian profiles to save CPU time. In addition, the NLTE models include 223,563 
secondary NLTE lines as well 219,774 LTE background lines, these hues were selected using 
the same criterion that we use for the LTE models. The calculation was performed with 
a variable resolution wavelength grid with 204,434 points, including extra points that are 
inserted to resolve NLTE lines. 

We show a comparison between the NLTE model and Vega observational data in 
Fig. 9. In general, the fit is good, no attempts have been made to fine-tune some of the 
abundances (e.g., to reduce the too strong absorption feature just shortward of I6OOA). A 
more detailed analysis of the Vega data will be presented elsewhere (Aufdenberg et al, in 
preparation) . 

We compare the NLTE and LTE structures for a Vega model with T^s = 9500 K in 
Fig. 11 and the departure coefficients of the ground states of some selected ions are shown 
in Fig. 10. The departures from LTE are larger for the Vega model than for the solar model 
shown above. The changes in the temperature structure are not very large in the deeper 
layers, but can reach nearly 1000 K in the outermost, optically thin, zones of the model. In 
the line forming regions, NLTE effects are generally small so that the spectra of the NLTE 
and the LTE models are not very different. Since the departures from LTE are larger for 
Vega than for the solar model, the validity of the assumption of LTE is weaker for Vega 
than for the solar model. This is the main reason why the LTE models that we present in 
this paper essentially stop at T^s — 10000 K, for higher effective temperatures NLTE effects 
will rapidly become important. 



4. Summary and Conclusions 

In this paper we presented a new set of models for main sequence stars for effective 
temperatures from 3000 K to 10000 K. We include mostly LTE models, however, NLTE 
effects become progressively important for effective temperatures larger than about 7000 K 
Therefore, we have included self-consistent NLTE models calculated with a subset of the 
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NLTE species that are available in PHOENIX. 

We find that our models agree well with the Kurucz 92 grid for solar type stars 
(5000 < Tefj < 7000 K) and that our LTE models for stars with effective temperatures 
up to about 10000 K also agree reasonably well with the Kurucz models. For lower 
temperatures we find some differences between the model grids. The reason for this is that 

our calculations arc designed to include the atmospheres of very low mass stars and brown 
dwarfs, so that they include a more detailed molecular EOS and a better set of molecular 
opacity sources than the Kurucz 92 models. For higher effective temperatures, NLTE effects 
become important. Our results for Vega show that NLTE has an effect on the structure 
of these models. Therefore, detailed NLTE models are required for effective temperatures 
higher than about 7000 K. Furthermore, Aufdenberg et al. (1998a,b) have shown that 
the combined effects of line blanketing and spherical geometry are essential ingredients in 
model atmosphere models for effective temperatures ^ 18000 K. We will present spherical 
NLTE models for larger effective temperatures in a subsequent paper (Aufdenberg et al, in 
preparation) . 

It appears that the assumption of LTE is a useful approximation only in a very narrow 
range of effective temperatures near that of solar type stars. For cooler and even more so, 
for hotter stars, NLTE effects are important and should be included in model atmosphere 
calculations. 
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Table 1: Complete listing of PHOENIX (version 8) NLTE species. The tables entries are of 
the form N/L, where N is the number of NLTE levels and L the number of primary NLTE 
lines for each model atom. The calculations presented in this paper include all species listed 
in the tabic with the exceptions of Li, Ti, and Co. PETER: you say you use Li and Ti in 
the sun and Ti in Vega. 





I 


II 


III 


IV 


V 


VI 


H 


30/435 












He 


11/14 


15/105 










Li 


57/333 












C 


228/1387 


85/336 


79/365 


35/171 






N 


252/2313 


152/1110 


87/266 


80/388 


39/206 


15/23 





36/66 


171/1304 


137/766 


134/415 


97/452 


39/196 


Ne 


26/37 












Na 


3/2 












Mg 




18/37 










Si 




93/436 


155/1027 








S 




84/444 


41/170 








Ca 




87/455 










Ti 


395/5279 


204/2399 










Fe 


494/6903 


617/13675 


566/9721 








Co 


316/4428 


255/2725 


213/2248 








total: 




5346/60637 
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Fig. 1. — Comparison between PHOENIX (full line) and ATLAS9 (dotted line) model spectra 
in the UV for T^s = 10000 K, log(c/) = 4.0 and solar abundances. The PHOENIX spectrum 
has been smoothed with a Gaussian kernel to lOA resolution. 
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Fig. 2. — Comparison between PHOENIX (full line) and ATLAS9 (dotted line) model spectra 
in the UV for Tcs = 10000 K, log{g) = 4.0 and solar abundances. The PHOENIX spectrum 
has been smoothed with a Gaussian kernel to 20A resolution. We also show the PHOENIX 
spectrum at its original resolution scaled by a factor of 0.5 for comparison. 

Fig. 3. — Comparison between PHOENIX (full hue) and ATLAS9 (dotted hue) model spectra 
for Teff — 4000 K, log{g) = 4.0 and solar abundances. The PHOENIX spectrum was convolved 
with a Gaussian kernel of 20A width for wavelengths smaller than 1 //m and with a Gaussian 
kernel of 50A width for wavelengths larger than 1 /xm. 

Fig. 4.— Comparison between the NLTE and LTE structure of PHOENIX solar NLTE models. 
Tstd is the optical depth in the continuum at 1.2/1. 

Fig. 5. — Comparison between the Kitt Peak Solar Atlas spectrum (thick curve) and the 
solar NLTE model (thin curve), both reduced to a resolution of AA = lOA. No attempts 
were made to fine-tune the model. Terrestrial features (e.g., at about 0.75, 0.95 and 1.12 /im) 
were not removed from the observed solar spectrum. 

Fig. 6. — Comparison between the Kitt Peak Solar Atlas spectrum (thick curves) and the 
solar NLTE model (thin curves) for selected wavelength regions. Both spectra were smoothed 
to a resolution of AA = 0.05A. No attempts were made to fine-tune the model. 

Fig. 7. — Comparison between the NLTE and LTE structure of PHOENIX solar NLTE models 
as function of the height above Tstd = 1 

Fig. 8. — Departure coefficients of the ground states of selected ions as function of Tstd 
(the optical depth in the continuum at 1.2//) for the solar model with T^s — 5770 K and 
\og{g) = 4.44. 

Fig. 9. — Comparison between a PHOENIX Vega NLTE model (dotted/thick fine) and Vega 
absolute spectrophotometry and lUE observational data. The model parameters are given 
in the figure, the "Vega abundances" are taken from CasteUi & Kurucz (1994). The model 
spectrum has been convolved with a Gaussian kernel of GA half- width to facilitate comparison 
with the low resolution data. 
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Fig. 10. — Departure coefficients of the ground states of selected ions as function of Tstd 
(the optical depth in the continuum at 1.2/x) for a Vega model with T^s — 9500 K and 
\og{g) = 3.95. 



Fig. 11.— Comparison between the NLTE and LTE structure of PHOENIX Vega NLTE 
models. Tstd is the optical depth in the continuum at 1.2/x. 



